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In this work we study the quantum Toda lattice, developing the asymptotic Bethe ansatz method 
first used by Sutherland. Despite its known limitations we find, on comparing with Gutzwiller's 
exact method, that it works well in this particular problem and in fact becomes exact as fi grows 
large. We calculate ground state and excitation energies for finite sized lattices, identify excitations 
as phonons and solitons on the basis of their quantum numbers, and find their dispersions. These 
are similar to the classical dispersions for small ft, and remain similar all the way up to ft = 1, 
but then deviate substantially as we go further into the quantum regime. On comparing the sound 
velocities for various ft obtained thus with that predicted by conformal theory we conclude that 
the Bethe ansatz gives the energies per particle accurate to 0(1/A'^^). On that assumption we can 
find correlation functions. Thus the Bethe ansatz method can be used to yield much more than the 
thermodynamic properties which previous authors have calculated. 

PACS Numbers: 63.20.Ry, 63.20.+t, 05.30.-d 



I. INTRODUCTION 

The Toda lattice [§, introduced by M. Toda in 1967 
1^, is a chain of particles which interact with nearest 
neighbours with an exponential potential. The quantum 
mechanical Hamiltonian for a periodic Toda system of 
length N (i.e. n + N = n) is 
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(1) 



where the u„ are displacements from equilibrium sites. 
We have chosen appropriate units to remove h, m (the 
mass of the particle) and the length scale of the potential. 
The infinite system also has a linear term in the potential 
(to cancel the one in the exponential), but with periodic 
boundary conditions this vanishes. 77 is a measure of 
the anharmonicity and also of the scale of the quantum 
effects. The larger rj is, the more "classical" the system 
and the more harmonic the low-energy excitations. In 
the classical limit the parameter i] can be scaled out but 
in the quantum case this can only be done by introducing 
an ft 7^ 1 in the above equation. We shall occasionally 
write 



ft 



(2) 



so that the Hamiltonian can be rescalcd and rewritten as 
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The Toda lattice is interesting, classically and quan- 
tum mechanically, because it is the one example of a 
nonlinear lattice which can be solved exactly. Elemen- 
tary excitations are cnoidal waves, which are periodic 



waves analogous to the normal modes of a harmonic lat- 
tice, and solitons, which are travelling pulse- like solutions 
which retain their shape even after interaction with other 
excitations. The periodic system does not support soli- 
tons of the infinite-chain type, since these involve a net 
compression, but a cnoidal wave with large amplitude 
behaves very much like a soliton (Fig. |]). 

The classical periodic system was studied by Kac and 
van Moerbeke, and Date and Tanaka It is com- 

pletely solved 0: given any initial condition of the sys- 
tem its future time evolution can be written down ex- 
actly. In quantum mechanics, there have been several 
treatments based on various approximations and assump- 
tions. On the one hand, Gutzwiller Q has given an exact 
treatment of the 3 and 4 particle lattices, and his quan- 
tization algorithm is capable of generalization to larger 
N as well. His results were rederived in the r-matrix 
formalism by Sklyanin and by Pasquier and Gaudin 
0. The method is plausible and makes a transparent 
connection with the classical formulation of the prob- 
lem. On the other hand calculating with this algorithm 
is a formidable task. The method is summarized in §IV. 
On the other hand, some authors have used the (asymp- 
totic) Bethe ansatz to treat the problem. Sutherland ||] 
originally recovered the classical results (high r]) in the 
thermodynamic limit {N ~~> 00). Later authors |^ have 
remained in this thermodynamic limit, but have looked at 
arbitrary 77, and have calculated various thermodynamic 
functions. 

In this paper, we use the Bethe ansatz to look at the 
case of finite N, which in some ways is more illuminating 
when one tries to classify excitations as phonons or soli- 
tons. §p obtains the Toda lattice as a limiting case of the 



1/sinh model, an idea due to Sutherland, and §[H sets 



up the Bethe Ansatz equations for the latter model and 
performs the same limit to obtain equations describing 
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the Toda lattice. Though the asymptotic Bathe ansatz is 
in general inaccurate for finite we find in this model 
that it is much better than it has been given credit for, 
and in particular becomes exact not only for N ^ oo but 
also for 77 ^ with finite N. In §[V we demonstrate this 



by setting up exact equations using Gutzwiller's method 
and seeing what approximations are involved in going 
from these to the Bethe equations. The claim that 
the Bethe ansatz misses a fixed fraction of states does not 
stand scrutiny. One need only glance at the harmonic 
limit (§^) where every one of the states is accounted for 
accurately. 

Having done this, we examine the opposite, highly 



quantum limit in §VI, which makes clearer how the low 



lying phonon-like modes go over to solito n-lik e states as 
their occupation number is increased. §VII calculates 



dispersion relations for phonons and solitons, and com- 
pares the classical and quantum results. We find that the 
results are essentially the same, apart from the quanti- 
zation of energy levels, for 77 > 2 (?i < 1) roughly; but 
as one decreases 77 further the quantum results deviate 
more and more from the classical, though they remain 
qualitatively similar down to 77 ^ 0. In this regime we 
get phonon-like excitations whose energies cannot be de- 
rived from harmonic approximations (why we think of 
them as phonons is discussed in §VII), and soliton-like 
excitations which can be thought of as authentic exam- 
ples of the much discussed "quantum soliton" . 



fiVIII considers finite size effects and makes contact 



with conformal theory to find correlation functions. We 
offer evidence that the asymptotic Bethe ansatz, in this 
problem, gives the energy per particle accurately to order 
l/N"^, though on general grounds it is guaranteed only 
to give results accurate to order 1. Finally, we consider 
in how all this relates to the classical lattice, and the 
appendix gives, for completeness, a brief discussion of the 
other conserved quantities (Henon's integrals) and why 
they are conserved in the classical and quantum cases. 



II. SCALING THE 1/ sinh^ MODEL TO THE 
TODA MODEL 

Sutherland was the first to treat the quantum Toda 
lattice, as a limiting case of the 1 / sinh^ model, by pi- 
oneering the use of the asymptotic Bethe ansatz. He 
contented himself with recovering the classical results, 
and showed that the classical solitons are recovered by 
taking the classical limit of particle like excitations of 
the Bethe equations. He did not explore regimes other 
than the classical, thermodynamic limit. Later authors 
like Mertens ||] have directly treated the Toda lattice 
by Bethe's ansatz, using the phase shifts obtained from 
the Toda potential; but the validity of the Bethe ansatz 
(which involves summing over phase shifts a given parti- 
cle suffers in collisions with all other particles) is unclear 





FIG. 1. How a classical cnoidal wave continuously goes 
from nearly harmonic to solitonic with increasing amplitude. 
Above, the first "normal mode" which goes into a one-soliton 
state. Below, the second mode which goes into a two-soliton 
state. The same thing happens in our quantum description 
when we put more and more phonons into a particular phonon 
mode. (Not to scale) 



in a model where only nearest-neighbour interactions ap- 
pear. We therefore use Sutherland's approach and scale 
the 1 / sinh^ model. Our scaling procedure is somewhat 
more explicit and displays the fact that the limiting pro- 
cess leaves us with a one parameter model, the Toda 
lattice with a general coupling constant 77 ( see below), 
from which the classical, the harmonic, and the extreme 
anharmonic limits follow. 

Our starting point is the Hamiltonian 
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■i^ sinh [{xra - a;„)/2a] 



(4) 



Here a is a length scale giving the range of the poten- 
tial, (7 is a coupling constant and the particles are on a 
ring of length L (so that the density d = N/L). In the 
dilute limit when the particles are far apart, the sinh^ 
becomes an exponential; we achieve this limit by making 
the substitution 



Xn = nld + u„a 



(5) 



where Un are displacements from lattice sites spaced 1 /d 
apart. We let ad go to zero, and assume that the u's are 
bounded (that is, the wavefunction vanishes as u — > 00). 
Then we have for m < n and da <C 1, 



sinh 



2a 



= sinh 



= (1/4) cxp 



2ad 

n — m 



ad 



and the potential in the Hamiltonian becomes 



ad 



(6) 
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So on putting 



The total momentum is 



9 = 



(7) 
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and then allowing ad to go to zero, all terms in the in- 
teraction except the nearest-neighbour terms (i.e. n = 
m + I) are killed, and we finally arrive at the Hamilto- 
nian 



N 



N 



,-(lln+l-tin) 
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which is the Toda Hamiltonian, eq. (m 



III. SOLUTION BY THE ASYMPTOTIC BETHE 
ANSATZ 

The 1/sinh^ system, being integrable, is characterised 
by N commuting integrals of motion. If we suppose that 
the particles are moved far away from one another, they 
don't interact except during short-range collisions, and 
for the rest of the time they have well-defined momenta 
which can be taken to be the conserved quantities. Dur- 
ing two body collisions the most that can happen is an 
exchange of momenta, and one can show that n-particle 
collisions can be completely described in terms of succes- 
sive 2-particle collisions and their phase shifts, so that 
the momenta are reordered but not changed. The Bethe 
Ansatz is a sum of plane-wave product states, charac- 
terized by a set of single-particle momenta {pn} and an 
amplitude for each plane wave state which features a dif- 
ferent permutation of these momenta. All that is required 
for calculation is the two-body phase shift 9{p — p') for 
two particles with momenta p and p' |^ ; then equations 
for the Pn can be written down and solved. These equa- 
tions are 
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Here /„ are integers for odd N, or half-odd-integers for 
even N, no two of which are equal. The energy (eigen- 
value of Hs) is then given by X^-P"^, so that the energy 
of the corresponding Toda problem would be o,^Pn^- 

This solution is derived in the limit when the particles 
are far apart, weakly interacting and in approximately 
plane wave states, so naively one would not expect the 
results to hold at higher densities. It is known, however, 
that this "asymptotic Bethe ansatz" holds at all densities 
in the limit N oo (the thermodynamic limit) provided 
the virial expansion has no singularities as a function of 
d PI . For this particular problem, it turns out that 
the solution is also exact for arbitrary N in the limit 
r] —> 0. Otherwise, though not exact, it is often a very 
good approximation. 



Ptot 



(9) 



Owing to the Galilean Invariance of the 1 / sinh^ model, 
a state with zero total momentum can always be boosted 
to have a total momentum ptot by adding an appropriate 
integer to all the quantum numbers, at a net energy cost 
plot I N independent of the coupHng constant. We can 
take the expression for the two-body phase shift in the 
1/ sinh^ system from Sutherland: 

e{p) = 2[argr(l + S + ipa) - argr(l + ipa)] (10) 



where S{S + 1) = 2ga'^ = ^e^ 



1^ — ■|e^/'"'. Since we are taking 
the dilute limit, 5* — > oo for any value of 7y, and we can 
write 



In the limit S* — > oo, the phase shift ( p^ becomes 
d{p) = 2palogS - 2 Im logr(l + ipa) 



(11) 



(12) 



[we can show this by using Stirling's expansion for large 
S in the first gamma function in ([lO|)]. We substitute for 
S from ([Tl|), put the resulting phase shift into the Bethe 
equations (|), noting that J2,n(Pn - Pm) = Npn - ptot 
where ptot is given by (^), and rearrange [the p„ on the 
left of (H) cancels with a term from the phase shift, leav- 
ing only 0{1/ L) and smaller terms]. Defining dimension- 
less "momenta" by fc„ = p„a, dividing out the common 
d and taking ad — > we end up with the equations to be 
solved: 



1 ^ 

— ^ Im logr[l + i{kn - A:™)] 



where for convenience we have written 
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logh). 



(13) 



(14) 



Note that the total momentum of the system, fctot = 
{ad){2-K / N) Ij, goes to zero as ad goes to zero, so that 
we are working in a zero momentum frame. This is a con- 
sequence of the length of the underlying 1 / sinh^ model 
going to infinity (on the scale of the range of the poten- 
tial) , the momentum being inversely proportional to the 
system length. However, our simultaneous scaling up of 
the interaction by an exponential factor (^ ensures that 
the individual particle momenta remain finite. Thus we 
have gone from a 1/ sinh^ "gas" with particles described 
by actual position coordinates, to a lattice with particle 
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positions described as displacements from lattice sites, 
and no net momentum, which is what we wanted. The 
energy of this Toda problem is J2^n- Since the problem 
continues to be Galilean invariant, a finite momentum 
fctot can always be introduced into the above equations 
by adding aktot/N to the right hand side, at a total en- 
ergy cost of k^^^/N. This fctot need not be quantized, 
since as the length of the underlying 1 / sinh^ model ex- 
pands the quanta of momentum become infinitesimal. 

The /„ in (||) and ( |l^ ) are the quantum numbers of 
the system, and uniquely specify the state of the sys- 
tem. The momenta kn are ordered in the same way as /„ 



[despite the apparently opposite sign for 77 > 2 in (13)] 
and we assume that the order is ascending in n. In the 
ground state the /„ are successive integers (or half inte- 
gers), generally taken to be centred about zero [though 
it does not matter here, since one subtracts their average 
value in Eq. (|l^)] and in the excited states one or more 
of them are increased by various integer values, always 
making sure no two of them have the same value. 

Although we took the dilute limit in arriving at these 
equations, the Toda Hamiltonian (^ which they describe 
contains no reference to the lattice constant, and there- 
fore they are valid at all densities, or at least at all den- 
sities sufficiently low that the particles do not cross each 
other. (The wavefunction will give the typical "spread" 
in Un and we must assume, for physical reasons, that 
the inter-particle separation is much larger than this). 
Mertens' treatment [^, if followed through, gives the 
same equations as the above but with an extra term on 
the right hand side equal to (27rad/iV) ^ /„ (which is the 
above ktot- he does not consider a limiting case of the 
1/ sinh^ model and does not take d — > 0). This term has 
no significance and, in particular, must not be confused 
with the phonon or soliton momenta (§V1]). In fact, it 
may be subtracted out, since it is independent of n, to 
recover our equations. We prefer this, the rest frame, 
because it is the frame in which one normally discusses 
phonons and also because it is convenient in making con- 
tact with Gutzwiller's work. 

Since one can add a constant quantity to the /„ 
without effect on the equations, they contain some 
redundancy — N — 1 quantum numbers are enough to 
characterize the system. We could define new quantum 
numbers by 



N-n+l 



N-n 



1, n=l,2,...,N-l (15) tion as * = C*-^! 



to spurious solutions where the ordering of the fc's is not 
the same as that of the Ps. 

Alternatively, one could pass to the thermodynamic 
limit and write down integral equations from which var- 
ious thermodynamic quantities could be calculated, as 
in Yang and Yang's treatment of the ^-function Bose 
gas. This has been done by Mertens, and by Hader and 
Mertens We define {N/2Tr)^{k)dk as the number of 
fc's between k and k + dk. Then (|l^) yields the integral 
equation for the density of the fc's in the ground state 
which is, in agreement with Mertens, 

^{k)^-2a+- [ (,{k') Re tpil + i{k - k'))dk' . (16) 

TT J-B 

(ip is the digamma function). For reasons given in the 
next section, Matsuyama []l2| gets the same equation for 
the distribution of the zeros of Hill's determinant in the 
Gutzwiller method (but without the inhomogeneous part 
since he takes ?i = 1, or 77 = 2.) 



IV. COMPARISON WITH GUTZWILLER'S 
FORMULATION 

The Bethe equations for the Toda lattice can also be 
derived from Gutzwiller's solution of the problem, if some 
approximations are made. This helps clarify what the k's 
mean in the non-dilute limit — in particular, their corre- 
spondence with the classical variables, and also tells us 
when our approximations are valid. We briefly describe 
Gutzwiller's method and the resulting quantization con- 
ditions. 

Gutzwiller, following the classical ideas of Kac and 
van Moerbeke |3j , tries to write the wavefunction of the 
N body lattice as a series involving the wavefunctions 
of the — 1 body open lattice obtained by removing 
one particle. Suppose these A^ — 1 body wavefunctions 
are vE'kik2...k;jv-i ; the indices k correspond to the classi- 
cal variables (the eigenvalues of the truncated A^ — 1 
dimensional Lax matrix). For the open chain they are 
purely imaginary but when using them as a basis in the 
closed chain Gutzwiller shows that one must extend them 
to have a real integer part; in other words, Ki = ipi + ki, 
where ki is an integer. One aims to find the spectrum 
of the p's. It turns out that if one writes the wavefunc- 



where the sum 



so that the Vn may take any integer value from upwards. 
(These are the number of "holes" between successive in- 
tegers /„, starting from the right.) These are, in the 
harmonic limit, the phonon occupation numbers (§0). 

Equations ( p^ can be solved numerically, for instance 
by the Newton-Raphson method, for moderate values of 
A^ without much difficulty if one has a good starting 
guess. If not, the numerical methods tend to converge 



is over the integers ki, one can get a solution of the form 



_i = (ki - K2)(k2 - K3) 



vided the coefficients r, s, 
relations 



■ r^^s^^tf,^ . . . pro- 
satisfy identical recursion 



-N 



(17) 



where D{k) is basically the characteristic polynomial of 
the Lax matrix (see the appendix): 
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d{k) = h'^K 



Eh 



N-2^N-2 

til 



... + {~if-^AN^inf, + {^ifAN, (18) 

and h is defined in (Q). Suppose that its zeroes are 
iei, i£2, ■ • ■ , if-N', then D{k) can also be written as 



D{k) = Y].{hK - iCn). 



(19) 



The same recursion relations are derived by Sklyanin, 
and by Pasquier and Gaudin, from different points of 
view They have two independent solutions, differ- 

ing in their behaviour at +00 and —00. Gutzwiller sets 



.(2) 



(-l)V" 



(20) 



where r' and r" are coefficients to be matched later when 
"joining" the two solutions, and r'^ and r" are two new 
variables which (it turns out) are complex conjugate. 
They have solutions 
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(22) 



-00, and the 



The former approaches a constant as k 
latter approaches a constant as k —00. 

If one tries to join these solutions, one gets the consis- 
tency condition 
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Here the -I- signs are for even N, the — signs for odd N. 
This determinant has N purely imaginary zeroes, which 
we call ipi, ip2, . . . , ipN (in ascending order). (If N is odd 
and all odd integrals ^3,^5,... vanish — this happens, 



for instance, in the ground state — then there are only 
— 1 zeroes but in that case k = also satisfies the 
quantization conditions below, so we include it among 
the p's.) It is clear that in addition to these, ip„ -I- 
where / is an arbitrary integer, are also zeroes of the 
determinant. 

The determinant is part of what we need to find the 
spectrum of p, but it is not enough since we don't know 
the constants of motion in D{k). We need more quanti- 
zation conditions; to supply these Gutzwiller defines an 
angle = (1/2) arg(r'/r") = arg(r') since r' and r" are 
complex conjugate. If one normalizes the solutions by 



ip 



(2) 

r- 
ip 



1, one finds 



= arg(r') = Im log 



ip 



(24) 

Then is a monotonically increasing function of p. Ab- 
breviating (pipn) as 0„, Gutzwiller's quantization condi- 
tion reads 

<Pi — 4'2 = ■ ■ ■ = (pN (modulo tt). (25) 
In addition he assumes that 

</.l+02 + -..+0Ar=O. (26) 

If both of these conditions are satisfied, the allowed values 
of (j)n are very limited; they can only be of the form /n7r-f 
mn/N, where m is the same integer for all n and /„ is 
an arbitrary integer, different for different n. But 0„ is 
an increasing function of p„, hence if the p„ are ordered, 
we must have the /„ also in increasing order. Then, from 
d^), we get 



■nin + "miT = 



(27) 



which yields m = — ^/„. So we have, finally, expres- 
sions for Gutzwiller's phase angles: 



N 



t- iNp„ 



argr = arg 



n™r(i + *[p„ -£,„]) 



-aNpn -f- ^ Im logr(l -t- i[pn - e™]) 



Im logr^ 



ip„ 



(28) 



or 



apn = 



TT 



In- 



N 



iogr(i + z[p„ 



N 

m 

Im logr^ 
iV 



(29) 



5 



{a = ilog(?7/2) = — log?i.) These, then, are the exact 
Gutzwiller equations which can be combined with ( |2^ ) 
to calculate the p„ and e„; once the latter are known, 
all the conserved quantities can be found. The /„ in this 
equation are the quantum numbers of the system, and 
are the same as the /„ in the earlier, very similar Bethc 
ansatz equations (13) — to which these equations in fact 



reduce provided (1) the last term can be ignored and (2) 
Pn is very close to e„ for all n. These things can happen 
under two circumstances. 

There is an argument in showing that the pn 
should approach e„ as ^ oo (and one knows on gen- 
eral grounds that the asymptotic Bethe ansatz is correct 
in this limit). This also happens as h —^ oo, for finite 
N. We can understand the latter fact intuitively as fol- 
lows: As ?i ^ oo, the polynomials D{k) tend to infinity. 
(This is not obvious — for example, they don't vanish as 
h ^ — but it will be demonstrated in §V]). Then they 
will be small only in a small region close to their zeroes, 
so the matrix of which ( |2^ ) is the determinant tends to 
the unit matrix except when p„ lie in some small regions 
surrounding e„. Thus the determinant can only vanish 
when the p's approach the e's; otherwise it is close to 
unity. For the same reason, (|l]) tends to unity (its ze- 
roes will be close to ipn + I ~ 0, I > 1, and for k — ien 
all the Cs in the denominators will be very large). Then 
the last term in ( |29| ) will vanish, and all the p's can be 
substituted with e's, and we recover exactly the Bethe 
ansatz equations. 

Thus the Bethe ansatz is actually more accurate in 
the quantum limit than in the classical limit. Indeed, 
even for Ti = 1 and A^ — 4-6, the agreement with Mat- 
suyama's exact diagonalization results |l3) is excellent 
(one gets exactly his answers, to his reported accuracy) 
— and this looks neither like a thermodynamic limit nor 
like an extreme quantum limit. 

One might imagine that the Bethe ansatz equations 
could be improved by subtracting the term Im log r^^ , 
but it turns out that this term is always small compared 
to the others and does not greatly improve the results, 
while it is computationally expensive to include; there- 
fore we ignore it in all cases. 

Finally, we observe that equations (|2^) do not remain 
the same if the p's and e's are increased by a constant 
quantity, because of the last term which does not ap- 
pear in the Bethe ansatz equations. We cannot therefore 
transform these easily to a non-zero-momentum frame. 



V. THE HARMONIC LIMIT (HIGH tj) 

For large 77 (the classical limit) the lattice is harmonic, 
at least for sufficiently small quantum numbers. The 
larger 77 is, the larger the energies and the quantum num- 
bers required for anharmonicity to show up. Treating 



this case makes clear the mapping between the phononic 
quantum numbers and the /„. 

First, the exact solution. There are A^ — 1 normal 
modes in the system, characterized by "phonon mo- 
menta" or wavenumbers qn = 2Tm/N, where n — 1 
. . .N — I. In our notation the coefficient of the terms 
in (^ is ri/2. Then the frequency cj„ of the nth mode is 



UJ„ 



2rysin(|) 



(30) 



An arbitrary state of the system is then characterised 
by a set of nonnegative integers { I'n } (phonon occupa- 
tion numbers). The energy of such a state is 



E = Nr]+ ^(j^„-f |)w„ 

N~l 

= Nrj + 2y/2^ ('^'^ + sin j 



(31) 



The first term arises from the constant term in the Taylor 
expansion of the exponential potential. For the ground 
state, we set = and find 



(32) 



which for large A^ has an expansion 



E = Nri+^/2^ [ -N 



6N 



+ 0{1/N'^)]. (33) 



Now we compare this expression with the result of solv- 
ing Eqs (|l^) numerically for 10 particles, and for various 
"large" f] : 10, 100, and higher. The ground state is 
when the !„ are contiguous, with no "holes"; the energy 
calculated from ( ^ij) is in good agreement with the value 
obtained from the Bethe ansatz. Numerical calculations 
show that the i/„ which describe a phononic state are ex- 
actly the numbers defined in (^5|). In other words, the 
number of phonons in a mode n is given by the number 
of holes between /jv-n and In-^u+i- 

This prescription accounts for all the states of the har- 
monic lattice, and the quantitative agreement is very 
close for low phonon numbers (the higher 77 is, the higher 
the allowed phonon numbers before anharmonic effects 
start showing up). Fig. |^ gives the dispersion curve for 
single phonons; only for low rj does it differ from the 
harmonic-lattice curve. Calculations show that the ener- 
gies of phonons are additive (provided there are not too 
many of them) so multi-phonon states are also accurately 
described. 



VI. THE STRONGLY-QUANTUM 
ANHARMONIC LIMIT (r/ ^ 0) 

In the large-77 case, increasing occupation numbers 
will bring out anharmonic corrections in the energy, and 
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modes with very high occupation numbers will resemble 
solitons. In §VII we demonstrate this with calculations, 



where 



but if 77 is not large anharmonicity shows up even in low- 
lying modes. 

Having looked at the harmonic limit in the last section, 
we now look at the opposite limit of the lattice, 77 — > 0; 
in this case it turns out that the phase shift simplifies 
greatly, and we can in fact solve equations ( ]l3| ) for fc„ — 
an uncommon phenomenon in Bethe ansatz calculations. 

Equation ( p^ ) for the phase shift is 

9{k) ^2klogS-2 Im \ogT{l + ik), 

and as — s- 0, fc also becomes small. In this limit the 
term involving the gamma function becomes 2jk, where 
7 = 0.577215... is Euler's constant. A quick way to derive 
this result is to assume S" is a large integer in (^0|) and to 
expand the first gamma function as a product {S+ik){S— 
1 + ik) . . . (1 + jfc)r(l + ifc), and if fc <C 1, the argument of 
this is k/ S + k/{S — I) + . . . + k + a piece which cancels 
the second term in (p^. As S* ^ 00, using the definition 
7 = lim„^oo 1 + (1/2) + (1/3) + . . . + (1/n) - logri, the 
phase shift becomes 



0(fc) = 2fc(log5 + 7) 



(34) 



(this is actually correct to quadratic order in fc), which 
when substituted in (||) yields 



kn 



27rd^^ + ^(7 + logs') J2 (^" - 



N 



N 



^In + 77 (7 + log S) iNk„ - fctot) 



and on substituting for fctot from 
find 



and rearranging, we 



Mil) 
N J 



(35) 



(Note that for very small 77, a will be large and nega- 
tive, so the negative sign above is deceptive; the fc's are 
ordered in the same way as the J's.) Equations ( ^5|) thus 
give fc„ for any excited state specified by any integers /„, 
and the energy is J2 ^ before. Note that the sys- 
tem now looks like a free Fermi gas or a hard-sphere gas, 
which indeed is the underlying model behind the asymp- 
totic Bethe ansatz (we derived our results as a limiting 
case of a gas of particles interacting by a 1/ sinh^ poten- 
tial). There is a continuous transition from this system 
to the classical Toda lattice as rj is increased. As we 
show below, even in this limit the excitations retain their 
qualitative features. 

In the ground state, the /„ are contiguous and may be 
taken to be 1, 2, . . . , N. Then a simple calculation gives 
the ground state energy as 



A 



iV2(a + 7)2 



(37) 



Now we consider excitations in which the last I /„ are 
excited by an amount m — we insert m holes between 
In-i and In~-i+i, or in the phonon language, we add 
m phonons in the Zth normal mode. J„ are now 1, 2, 3, 
...,N-l,N-l + m + l,N-l + m + 2, ...,N + m + l. 
Again, one can calculate the excitation energy; it is 



We consider several cases: 



E-Eo^ Aim - P) ( m + ^ 



(38) 



1. m small, arbitrary I 

In this case, we get approximately 

E-Eo= A{Nl - f)m. 



(39) 



This looks very much like a phonon dispersion; it 
rises from zero to a maximum at the zone boundary, 
where its slope dies off. It is linear in the number of 
"quanta" m, and for the lower-energy modes (lower 
I) it is also linear in mode number or wavcnumbcr 
(i.e. the 2nd mode has twice the energy of the first 
mode, and so on). 

Moreover, for phonons we know that the zero point 
energy in each mode is half the energy of one 
phonon; we can therefore sum half the above ex- 
pression over Z, for m = 1, and S66, RS 0j check, 
whether we recover the zero point energy ( ^6[ ) . And 
indeed, we do get 

N-l 

^ \A{Nl ~ f) = ^,AN{N' - 1) 

i-i 

in agreement with (|36|). 

The excitations are non-interacting — if we have 
several such excitations in different modes their 
combined energy is the sum of their individual ener- 
gies, if there are not too many of them. These hole 
excitations are thus quite analogous to phonons, 
though they cannot be derived by approximating 
the lattice to a harmonic lattice. 

2. / = 1, m large 

These are the excitations which one would expect 
to be soliton-like. In this limit, we get 



(40) 



E^Eo^AiN-l)[m+ — 



Eo = AN{N^ - 1)/12 « AN^ /12 



(36) 



For large m the energy is thus quadratic in m. 
This energy, however, is measured in the zero- 
momentum frame which is not the frame in which 
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one normally discusses solitons. The question of 
what is the correct frame is discussed in the next 
section, where dispersion relations are derived. 

I small, m large 

From ( |38| ) we note that if I is small, the excitation 
energy is proportional to I. For instance, the en- 
ergy for ? = 2 is twice that for 1 = 1. It is tempting 
to suppose that this is a two soliton state, since the 
energies of solitons are additive provided that they 
are few in number and hence well separated "most 
of the time" . In that case there would be a contin- 
uous transition between a phononic excitation of 
the second normal mode and the two soliton state, 
just as there is between the excitation of the first 
normal mode and the one soliton state. (Cf. Fig. 
S and §0) 



If the last two integers are excited by different amounts, 
one would presumably have two solitons with different 
energies. Here, too, the total excitation energy is the 
sum of the individual energies. Carrying this picture fur- 
ther, an — 1 soliton state (with all solitons having equal 
energies — / = iV — 1, m large) has all the particles except 
one moving in one direction like hard spheres, and is re- 
lated by a Galilean transformation to a 1 soliton state. 
An N soliton state (with all solitons identical) is simply 
a uniform translation of the lattice as a whole. One can- 
not put more than N solitons in an A^-particle lattice. 
The last few sentences are speculative, but they indicate 
the possibility of writing an arbitrary excited state as a 
kind of nonlinear superposition of solitons. (To make this 
more convincing, read cnoidal waves for solitons). Much 
the same thing is done in the classical periodic system 



VII. 



DISPERSION RELATIONS FOR PHONONS 
AND SOLITONS 



We now find the dispersion relations for phonons and 
solitons. First, however, we clarify the meaning of the 
momentum of these excitations. 

As remarked earlier, the fact that we take the dilute 
limit gives us a zero total momentum. Mertens' treat- 
ment, on the other hand, yields a finite momentum ^ fc„ 
proportional to ^ /„ and to the density d. This momen- 
tum is not a physically relevant quantity. It is not the 
momentum of a phonon (though it is proportional to it), 
since it depends on r; while the phonon momentum is a 
purely geometrical quantity depending only on the sys- 
tem size and lattice spacing. Nor is it the momentum of 
a soliton (it is not even proportional) since the soliton 
momentum doesn't depend on the lattice spacing. 

The phonon momentum q is the wavenumber of an 
oscillatory excitation. For an N particle lattice q has N 




FIG. 2. Phonon energies (for the Hamiltonian (j^ in units 
of h), plotted against wavenumber q for various rj. The solid 
line is the harmonic lattice curve and the curves for all "large" 
values of ri lie on top of it. The dotted line is for rj = 2, 
the dashed line rj = 0.1, the dotdash line rj = 0.01. The 
range of q is [— 7r,7r]. Energies are in units of (using the 
Hamiltonian (^); or with the Hamiltonian (^) energies are in 
units of h). Here and in later graphs, units are chosen to get 
an r;- independent curve in the large rj limit. 



equally spaced values generally taken to lie between — tt 
and TT (the first Brillouin zone) in units of the inverse 
lattice spacing. The soliton momentum is a little trickier 
to define in the quantum case. We discuss it below. 

First consider the small rj limit. We consider a single 
phonon, occupying normal mode n; its excitation energy, 
from (II) with / = n and m = 1, is E-Eq = A{Nn-n'^) 
and its wavenumber q, in units of inverse lattice spacing, 
is 2TTn/N (modulo 2tt; we can choose the value to lie 
between — tt and tt.) Note that Im = n in this case, 
if it was taken to be zero in the ground state, so q is 
proportional to this quantity. This gives w, the frequency 
(or the excitation energy of one phonon, since ?i = 1) in 
terms of q as 



UJ 



q 



2(7-ha)2 
and the phase velocity of sound is 



< g < 27r 



TT 



" 2(7-ha)2 
while the group velocity is 



1 



2tt 



2(7 + a)2 



(41) 



(42) 



(43) 



(in units of the lattice spacing). 

In the classical limit, of course, the phonons are what 
one would find from a harmonic approximation. For a 
mode with wavenumber q the energy is 



E-Eo = 2y2^sinig 



(44) 
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q (wavenumtaer) 



FIG. 3. The velocity of sound, dE/dq, plotted against q 
for various 77 for a 19 particle lattice. The solid line is the 
curve for the harmonic lattice, valid for large r). The crosses 
represent rj = 2, the circles r] = 0.1, the asterisks 77 = 0.01. 
Units are as in the previous graph for the phonon dispersion. 



which yields the phase velocity (in units of lattice spac- 
ing) 

, — sin i(7 

y = 2^ ^ 45 

q 

and the group velocity 

Vg = y^cosig. (46) 

The relations are different in the two cases, but have some 
similar features, and at intermediate values of r/ one ob- 
tains interpolations between these. Dividing the energies 
of excitation by -y/^ry one gets results independent of 77 
in the classical limit ?i — > or ry — s- cx). The results are 
plotted in Figs. || and || (for a nineteen particle lattice). 
One observes that for 77 > 2 the dispersion is more or less 
the classical harmonic-lattice dispersion, while it begins 
to deviate for ?/ < 2. This is further emphasized by Fig. 
^ which shows how the long wavelength sound velocity 
varies with a — — log?/. 

When we consider a soliton, we have to make clear 
what frame to view it in to obtain an appropriate mo- 
mentum. In the classical case it is usually viewed in 
the frame where "most" of the particles are at rest and 
only a localized excitation is moving. We would like to 
choose a frame in the quantum case such that the dis- 
persion agrees with the classical formula; in particular 
as the energy of the excitation increases it behaves more 
and more like a single hard sphere moving in a stationary 
background and the energy tends to fc^ (plus the ground 
state energy). 

We identify a soliton with a state where is greatly 
excited compared to all the other /c's. We can achieve 
the fc^ dispersion if we work in a frame where the fc's ex- 
cluding kN are (roughly speaking) centred around zero. 
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-6-4-2 2 4 6 

alpha 

FIG. 4. Variation of long-wavelength sound velocity, in the 
same units as in Figure H, as a function of a (= — log h = 
(1/2) log(r;/2)). Note that a — seems to divide the harmonic 
and quantum anharmonic regimes, i.e. the region where the 
harmonic approximation is valid for small excitations and the 
region where the zero point motion is so large that the har- 
monic approximation is not valid even in the ground state. 



In that case for large excitations kp^ ^ kn {n < N), 
the total momentum is very nearly kN, the total en- 
ergy is nearly fc^r^ « k^, and the quadratic dispersion 
is achieved. 

However, exactly how to define the frame is not clear. 
There are various possibilities — one could choose the av- 
erage of all kn except fci and kN to be zero (so that the 
/c's are not very much displaced from the ground state 
value); one could make the average of fc„ including ki 
but excepting kN zero; one could fix one of the fc's (say 
ki, kN/2 or kN-i) to its ground state value; and so on. 
These possibilities are plotted in Fig. ||, for 77 — 1000, and 
the dispersion for a classical cnoidal wave of wavelength 
N plotted for comparison, calculated from the formula 
for Tn = Un — Un-1 givcu in ||l[ (cf. §p^). Of the pos- 
sibilities listed the second (where the fc's excepting kN 
average to zero) seems the closest to the classical curve, 
but the agreement is imperfect and the "correct" frame 
would appear to be something close but slightly differ- 
ent. In plotting these curves we have used the Hamilto- 
nian (^, whose limit as ?i — > is the classical problem in 
the correct units. Fig. || shows the dispersion curves for 
77 = 1000. 

Fig. ^ shows the particular dispersion curve obtained 
by averaging fcn<7v to zero, for various 77. As in the case 
of the phonon curves, the soliton dispersions lie on top of 
each other for large 77 but begin peeling apart for 77 « 2; as 
rj is reduced further they move further and further away. 
Thus we find again that 77 = 2or?i=lisa boundary 
between classical and quantum regimes. For higher 77 the 
dispersions are essentially the classical ones apart from 
the discreteness of the energy levels. For lower 77 the 
results deviate significantly from the classical ones. All 
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FIG. 5. Dispersion curves for the classical and quantum 
cnoidal waves. The solid curve is the classical cnoidal wave 
or soliton, viewed in the appropriate frame. The dotted 
curve is the quantum cnoidal wave in a frame in which 
kn = 0. This lies closest to the classical curve among 
the cases considered. The dashed curve corresponds to the 
frame in which only the "inner" fc's are centred at zero, i.e. 
'^^_2 kn = (as in the ground state of the system). The 
dotdash curve corresponds to fc]v — 1 being fixed at its ground 
state value — implausible perhaps but included here for vari- 
ety. All curves are for = 10 and (in the quantum case) 
for ?7 = 1000. The Hamiltonian (|) is used; in other words 
energies from the Hamiltonian (llh are plotted in units of rj. 



the curves above have been calculated for a ten particle 
lattice. 

In the 77 ^ limit we have the fc„ given by (|3^) ; for the 
ground state state we take the /„ to be centred at zero 
(i.e. they range from -{N - l)/2 to {N - l)/2 for odd 
N, or from —N/2 to N/2 for even N); and for the soliton 
we excite Ipf by an amount m. Then ^ /„ — m. Clearly 
if we want the fc„ (for n < N) to be centred at zero, we 
must add to ( |35| ) a quantity to cancel the J2 In/N in the 
numerator, and instead subtract ^n=i -^n/i^ ~ !)• In 
this new frame, we have 



FIG. 6, 

1 "'rj 



kn — 



N-1 



Ni-f + a) 
^ " Nij + a) ' 

rr? + Nm + 



iV2(7 + a)2 



N{N ^\) 



(47) 
(48) 

(49) 



The energy formula is not very different from (|4C|). The 
details of this formula should not be taken very seriously 
since we are not clear about what the appropriate frame 
is in which to view the soliton. But the essential idea, 
that the energy is quadratic in the momentum at large 
energies, will remain. In this frame the energy is in fact 



The soliton dispersions, plotted in a frame in which 
kn = 0. For 7; > 2 all the curves lie on top of each 
other; they are shown by the solid line. The dashed line is for 
= 2 when they just start peeling apart. The dashdot line 
and the dotted line are for 77 = 0.1 and 77 = 0.01 respectively. 
The energies are in units of 7;, i.e. the Hamiltonian (^, in 
terms of 7i, is used. 



(apart from a constant piece) purely quadratic in the 
momentum — there is no linear term. This can be recon- 
ciled to our picture of the low-77 limit as a hard sphere 
gas, so that at any time the entire energy apart from the 
zero-point contribution comes from the kinetic energy of 
one particle, the other particles being at rest. 

Finally, if we wish to compare our system to the free 
Fermi gas which it resembles in one limit, we could look at 
the "particle-hole excitation spectrum" commonly plot- 
ted for such systems. To do this we start from the ground 
state, with contiguous /„; pick up one of these, say 
move it to /,'„ (where > /jv since all other states are 
occupied), and define the momentum of this "particle- 
hole excitation" as Q = 2tt{I^ — Im)/N. (This is basi- 
cally the total phonon momentum of such an excitation.) 
Then one gets a one-parameter range of energies for every 
Q, as shown in Fig. 0. The harmonic and rj limits 
look similar, qualitatively; the phonon or hole branch 
(the lower edge for Q < 2-k) is a sine curve in the for- 
mer case and a parabola in the latter, and the particle 
branch (the upper edge and the lower edge for Q > 2tt) is 
a straight line in the harmonic limit and a curve (which 
indicates nonlinearity) otherwise. The upper edge of the 
particle hole continuum has been identified with a "soli- 
ton" by Sutherland, and corresponds to promoting fcjv 
from the ground state configuration to one with a larger 
value, and is thus essentially identical to our picture ex- 
plained above. A study of the quantum numbers of the 
solitons and the phonons leads to a suggestive "phonon 
decomposition" of the soliton: we can view the soliton 
creation operator Aq schematically in terms of a phonon 



creation operator aj as 

aL 



(50) 
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FIG. 7. Particle-hole excitations. In between the bounding 
upper and lower curves lies a continuum of allowed energy 
values corresponding to each Q where Q is as defined in § VII 
for a single particle-hole pair. The upper graph corresponds 
to the harmonic limit, the lower graph to the r; — + limit. 
(The energy scales are different in the two graphs.) 



i.e. a particular kind of highly symmetric multi phonon 
state. 



VIII. CORRELATION FUNCTIONS, 
FINITE-SIZE EFFECTS AND CONFORMAL 
THEORY 

We now turn to the issue of correlation functions of 
the Toda lattice, making contact with the the theory of 
conformal invariance in this class of systems. Conformal 
invariance has given considerable insight into correlation 
functions of quantum many body models having critical 
behaviour, as typified by a vanishing of excitation ener- 
gies or power law correlations, and useful reviews of this 
fast growing field are to be found in and in [ p^ . 

Let us first note that the quantum Toda lattice in its 
ground state is not quite a lattice: the Bragg peaks are 
melted due to zero point motion. In the harmonic limit 
this is simple to see, since we can write the displacement 
in terms of the phonon creation operators and the phonon 
dispersion ujq — 2v\ su\{q/2)\ as 



exp(i9n) 



1 



(4 



(51) 



whereby «) = (l/iV)E(lH) ^ {l/^v)\og{N). The 
phonon velocity v — \/2/q in the harmonic limit of the 
Toda problem. The structure function at the first recip- 
rocal lattice vector G = {2tt/N) is 

{PG Pg) =^(e^2.n„ g-.2.n,„^ 



where we have used the Gaussian cumulant theorem 
(exp(a)) = exp(l/2(a^)) and the logarithmic integral 
(1/A^) ^(1— cos(gr))/a;q ~ (I/tt?;) log(r/ro). We thus see 
that the Toda lattice may be expected to have powerlaw 
correlations for all r/, since it has low energy excitations 
for all rj, namely the phonons. 

A characteristic of conformally invariant theories is the 
'central charge' c. One way of checking for conformal 
invariance is to compute corrections to the ground state 
energy for a finite sized system, which is expected to have 
a behaviour 



i?(L) = Leoo - ^ + 0(1/l2) 



(53) 



g-27r^{(u„-«„)^> ^ g-47r/i;log(|m-^n|) 



where v is the velocity of the low lying excitations, 
such that a tower of excited states exists with energy 
v2it/ L X integer. A glance at ( |3^ ) shows that in that limit 
of large rj we have c = 1, as indeed does the initial 1/ sinh'^ 
model. The case of c = 1 usually leads to exponents 
varying continuously with coupling constants, and hence 
( ^ ) is consistent with this possibility. In the present 
model, we must, however, first establish that the asymp- 
totic Bethe ansatz gives the correct energy to 0(1/A^) or 
0(1/L). This is not guaranteed a priori by any theoret- 
ical argument and must be checked for self consistency. 
(Incidentally, in the Toda lattice we are at a fixed den- 
sity so we will not distinguish between L and N .) The 
internal check performed is to compute the velocity at a 
fixed rj and to compute the energy for various N and to 
check against Eq. (|5^). 

First we note that in the extreme anharmonic limit 
equation ( |36| ) for the ground state in the low 77 limit does 
indeed give the same sound velocity as ( ^2|) or (|4^), so 
in the low 77 limit c = 1 exactly, as it is in the harmonic 
limit. 

We performed the calculation for rj = 2, 10, 100 (Table 
p. As in figures || and ^, we use the Hamiltonian ^ 
and units of h (equivalently, the Hamiltonian ([^) with 
units of y/2rj); in these units the sound velocity for the 
harmonic lattice is 1 exactly. From these results, we get 

vc^ 1.0764 ±0.0006 (r/ = 2) 
1.035 ± 0.003 (r/=10) 
1.01 ±0.03 (r/==100) 

On interpolating the 19-particle results of Fig. ^ for q=0, 
we get the estimates v = 1.08, 1.04, 1.01 for r? = 2, 10, 100 
respectively, with uncertainties in the second decimal 
place. Thus we get for the central charge 

c= 1.00 ±0.01 (ry = 2,10) 
1.00 ±0.03 (r/ = 100). 

The uncertainty in the cases 77 = 2, 10 arises mainly from 
the inaccuracy in our determination of v. The results 
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TABLE I. Ground state energy per particle as a function 



of system size 




N 




E/N 






r, = 2 


77 = 10 


77 = 100 


29 


1.675512397777 


2.890224040772 




33 


1.675665073759 


2.890370838890 




41 


1.675847391894 


2.890546128059 




49 


1.675947956073 


2.890642813964 




57 


1.676009234239 


2.890701728785 


7.713369630480 


65 


1.676049312911 


2.890740261674 


7.713407380691 


81 






7.713452036110 


97 






7.713476466168 


113 






7.713491273955 


129 






7.713500922044 



5E = 



'^tot 

N ' 



seem to indicate that c is equal to 1 at all values, and 
moreover it is reproduced correctly by the Bethe ansatz 
even at 77 = 100, which is well into the "classical" limit. 
It thus appears that the error in energy per particle goes, 
at worst, as the inverse cube of the number of particles. 
The error bars could be reduced by increasing the system 
size further. 

In the anharmonic limit, in fact, the series stops there 
[E/N has only a constant piece and a piece) while 

in the harmonic limit all odd powers, and 
so on are missing. One might conjecture that this is the 
case at all values of rj. For rj — 2 we took the ground 
state energies per particle for various N, subtracted Coo 
and the l/N"^ piece, and fitted the results to power series 
in 1/N starting at N~^. The result was a coefficient 
of 0.021 ± 0.006 for the N'^ term, and -1.1 ± 0.2 for 
the TV"'' term. Thus the coefficient of the term 
does seem to be very nearly zero. For 77 = 10 and 100 
the numbers we obtained didn't allow us to make such 
fits — the error bars turned out to be much larger than 
the values themselves. We conjecture that the coefficient 
of the term vanishes at all 77, but for high 77 the 

Bethe ansatz may not be accurate to this order in N and 
may be unable to reproduce this result. We are unable 
to make a statement about higher odd powers. 

Accepting that the Toda lattice is a c = 1 theory, we 
can establish the power law of the density correlator as 
in (^2|), without too much detailed calculation, on using 
the Galilean invariance of the model. The theory of con- 
formal invariance (see e.g. [|l5|) says that if we have an 
excitation that boosts the total momentum by ktot then 
the change in energy is 



SE = 2ttvx/N 

( ktot N 2 
X = (^) 

a — 2fj, 



(54) 
(55) 
(56) 



hence we find 



V 



(57) 



(58) 



where a is the exponent determining the decay of a pri- 
mary operator. However, Galilean invariance implies 
that 



Comparing with the harmonic limit result (p2|), we see 
that the primary operator may be identified with the 
density fluctuation pc and hence the result (|5^ ) is true 
at all 77 provided we substitute the appropriate value of 
v{r]) . A similar result is well known to be true for the 1 /r^ 
models for the density correlation function, but unlike 
in that case, there is a difficulty in defining a "bosonic" 
correlator, since we are always working at a fixed density, 
and hence the compressibility is zero. 



IX. COMPARISON WITH THE CLASSICAL KAC 
- VAN MOERBEKE FORMULATION 

To summarize the above, we now have a picture of how 
the fc's in the Bethe ansatz (or the p's in Gutzwiller's 
treatment) behave, in the ground state and in the ex- 
cited states. In the ground state the /'s and therefore 
the fc's are all closely spaced. In the excited states the 
separations between them widen. If there is a gap of 
m integers between In-u and Iff-n+i the gap between 
k^-n and k^-n+i widens and one has tti phonons in the 
77th normal mode. If the gap between the fc's becomes 
very large the excitation becomes solitonic. In particular, 
for 77 = 1 one has a one-soliton state; for 77 = 2, a two 
soliton state (with equal amplitudes); and so on. 

We now compare this description with the descrip- 
tion of the system in the classical variables of Kac 
and van Moerbeke ||,|l|. Briefly they use the variables 
/ii, /i2, . . . , fJ.N-1 which are the eigenvalues of a truncated 
Lax matrix obtained by striking off the flrst row and col- 
umn (i.e. removing the flrst particle from the problem). 
These /i's are the momenta of the particles in the re- 
maining open chain if the system is dilute. Kac and van 
Moerbeke show that these /x's are confined to the — 1 
closed intervals where the characteristic polynomial of 
the Lax matrix, |A/ — L\ is equal to or greater than 2 
in magnitude. The polynomial goes to ±00 for large A, 
while it oscillates in the middle; for the ground state it 
touches the lines A = ±2 in — 1 places so that the 
closed intervals referred to above are single points and 
all the /i's are stationary. For an excited state the poly- 
nomial crosses the lines A = ±2 so the closed intervals get 
a finite width and the fi^s oscillate inside these intervals 
as the system evolves. 

The analogues of the classical /x's are Gutzwiller's p's, 
or approximately Sutherland's fc's. Whereas there are 
iV — 1 T^'s each confined to a different interval in the 
classical picture, in Gutzwiller's picture each of the TV — 1 
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analogous variables has a spectrum of A'' values p„. On 
calculating the classical /z's in the ground state, as is 
done in , we find that their values lie almost exactly in 
between the quantum p (ie k) values. There is an analogy 
between the ^'s and the "gaps" in the fc-spectrum. In 
the ground state the gaps are minimum, the /u's fit into 
these gaps, and the /it's are stationary. In an excited 
state some or all of these gaps between the fc's widen, 
and the corresponding /i's are no longer stationary but 
oscillate in intervals of finite width. In particular a pure 
cnoidal wave corresponds to exactly one p, acquiring a 
width in which to oscillate, or exactly one gap among 
the In (hence the kn) widening. 

A single cnoidal wave has the formula j|] 

e"''" = 1 + {2Kvf (dn^ [2(n/A ± vt)K] - E/K) , (59) 

where r„ = u„ — Un~i, K and E are the complete elliptic 
integrals of the first and the second kind, A is the wave- 
length {— N for the first "normal mode" or one soliton, 
N/2 for the second normal mode, etc) and v is given by 



2Kv 



1 



E 

sn2(2if/A) ""^^^ 



-1/2 



(60) 



For low modulus k of the elliptic functions, this is like 
a sinusoidal wave, but as the modulus increases it be- 
comes sharply peaked locally and flat elsewhere (Fig. ||). 
As remarked in §VII, the dispersion calculated from this 
expression is close to the dispersion, in an appropriate 
reference frame, of the quantum cnoidal wave. 



X. CONCLUSIONS 

In conclusion, we have shown that the usefulness of the 
asymptotic Bethe ansatz in the quantum Toda problem is 
not confined to finding thermodynamic properties. The 
method gives results for energy per particle accurate to 
0(1/A^2), which is sufficient to calculate finite size effects 
and even correlation functions using conformal theory. 
The 0{1/N^) term seems to vanish in the exact solu- 
tion, though the Bethe ansatz solution probably doesn't 
reproduce this result. 

We have demonstrated that in fact the Bethe ansatz 
equations are a simplification of Gutzwiller's method and 
can be derived from them. The parameter governing 
the error can be taken to be the difference in p„ and 
e„ in § |v| . According to Matsuyama [|l2| this differ- 
ence falls exponentially with N , so that the error goes 
as e~^^^'^'^\ where /(r;) is some dimensionless number. 
We also show that the error vanishes as 77 becomes small, 
so that /(7y) ^ as 77 — > 0. 

Thus, we can treat finite sized systems, account for 
low-lying states (phonons) and higher excitations (soli- 
tons), and find their dispersions and velocities. Com- 
parison with conformal theory gives the "central charge" 



c = 1, which means that the coefficient of the l/TV^ term 
in the E/N expansion is essentially the sound velocity. 

We find that the properties of excitations are very sim- 
ilar to the classical properties for 77 > 2 (?i < 1), apart 
from the underlying discreteness of the energy levels. The 
quantization is then analogous to the quantization of a 
harmonic lattice. The soliton, which is an effect of large 
occupation of one mode, is no different from the classical 
object described by Toda; even its energy is effectively 
not quantized since the occupation number is so large. 

For small 77 (large Ti) things are different: the phonons 
no longer derive from a harmonic approximation, and 
the soliton dispersions no longer match the classical 
ones, though qualitatively the dispersion curves retain 
some similar features, both for solitons (high-amplitude 
cnoidal waves) and for phonons. For both excitations 
the dispersions depend on 77, and moreover the energy 
of a mode deviates rapidly from linearity with increas- 
ing occupation number 77, so that 77 need not be macro- 
scopic (at least for finite lattice size N) for the mode to 
become soliton-like — the soliton's energy is indeed quan- 
tized. Thus if the large 77 soliton is essentially the soliton 
of Toda's classical lattice, the corresponding small 77 ob- 
ject deserves to be called the quantum soliton. 



APPENDIX: HENON'S INTEGRALS, 
CLASSICAL AND QUANTUM 

In this appendix we discuss the integrability of the 
Toda lattice classically and quantum mechanically; while 
much of the discussion is not new it seems difficult to 
find it in one place elsewhere. Following Pasquier and 
Gaudin 0], who give a proof of quantum integrability, 
we show that their conserved quantities are the same as 
Henon's integrals — whose conservation is necessary for 
Gutzwiller's treatment to go through. 

The equations of motion for the classical lattice can be 
written in the Lax form 



dL 
'dt 



LM ~ ML. 



(Al) 



where 



L 



M = 



/ ^1 fli 

a\ hi a2 
0.2 h 

/ ai 
ai a2 
-a2 

V aN 



UN \ 



(A2) 



-a AT 



-OAr-i / 



(A3) 
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and 



-(9ji-<?3)/2 



Pj- 



(A4) 



From this one can show that the eigenvalues ol the 
Lax matrix L, or equivalently the coefhcients /„ ol the 
characteristic polynomial ol the Lax matrix are conserved 
quantities These are Henon's integrals, and are given 
by 



where 



il.'2 ifc. 



PtlPt2 ■■■P^k (-^Jl) i-Xj^) . . . i-Xj,) 



(A5) 



(A6) 



there are no repeated indices in the p's or the g's in a 
given term (zi, 12, . . . , ji, ji + 1, 32, j2 + 1, ... are aU 
different) , the total number of such indices in each term 
is m (i.e. fc + 2/ = m) and the sum is over all distinct 
terms satisfying these conditions (i.e. terms not differing 
merely in the order of factors). 

In quantum mechanics, the coefficients of the Lax ma- 
trix are the same, and have no ordering problems, but 
now the equations of motion (Al) are no longer valid 
(each term in the matrix product has to be ordered) 
and the proof that the coefficients are conserved fails. 
Gutzwiller [| assumes that they are conserved nonethe- 
less (he only takes the cases = 3,4 where it can be 
verified easily). Their conservation can be shown as a 
consequence of the work of Pasquier and Gaudin Q , who 
prove that the coefficients of u in the trace of the "mon- 
odromy matrix" Tjsi are in involution, where 



T/v(m) = L1L2 



Ln{u) 







(A7) 
(A8) 



The definitions hold in both the classical and the quan- 
tum cases. Classically their conservation follows from the 
classical equations of motion 



where A/„ 



u 





(A9) 
(AlO) 



Quantum mechanically these satisfy the Yang-Baxter 
equations: we may rewrite L„ — > Ln,g{u) = [u — p„)(l -|- 
o^g)!'^ ~ exp(— 5,1) + exp(g„) and show that the 
monodromy matrix Tn(u) Tg{u) satisfies the Yang- 
Baxter condition Tg{u)Tg'{v)Rg_g'{u — v) = Rg^g>{u — 



v)Tg>{u)Tg(v) with Rg 



a{u^v) + b(u — v)(7'g.cf^>. Tak- 



ing a trace over the auxiliary spaces ag , ag' the integra- 
bility is established. We now show that these coefficients 



are in fact Henon's integrals. Consider a polynomial in 
w, Fn{u), defined by 

...(-A,J(-X,g...(-A„) (All) 

where the indices satisfy the same restrictions as in the 
definition of Henon's integrals. It is easily seen that 



N 



,N~ 



(A12) 



We can show by induction that this polynomial is the 
trace of Tn{u). Defining 



F^{u) — All the terms in Fn{u) 

which do not include a factor e'™ , 
F^{u)^FMiu)~F;j{u) 

— All the terms in -Fjv(u) 

which include a factor e'" , 



we claim that 



T, 



z'^-F'^^^{u) 



N 



e-'i-+-F'^^M F'J,{u) 



(A13) 



(AM) 



(A15) 



The claim is easily verified for TV = 1, 2, etc. Suppose it 
is true for N\ then 



Tn+i — TnLn+i 



( {u - Pn+l)F'j^{u) 



29" + iFjV(u) \ 



{u-pN+i)e-'i-+-F'J,^,{u) F'J,^,{u) 
\ ~e-'i-^-F'^(u) 



which one can check is the same as 

F'^^M e'i«^-F'^{u) 



(A16) 



T. 



N+l 



e-'^-^-F^^^iu) F^^,{u) 



(A17) 



Thus our claim is true for all N, and in particular the 
trace of Tn is Fn{u). 
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